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Abstract 

Using a continuum description, we account for fluctuations in the ionic solvent surrounding a 
Gaussian, charged chain and derive an effective short-ranged potential between the charges on the 
chain. This potential is repulsive at short separations and attractive at longer distances. The 
chemical potential can be derived from this potential. When the chemical potential is positive, 
it leads to a melt-like state. For a vanishingly low concentration of segments, this state exhibits 
scaling behavior for long chains. The Flory exponent characterizing the radius of gyration for long 
chains is calculated to be approximately 0.63, close to the classical value obtained for second order 
phase transitions. For short chains, the radius of gyration varies linearly with N, the chain length, 
and is sensitive to the parameters in the interaction potential. The linear dependence on the 
chain length N indicates a stiff behavior. The chemical potential associated with this interaction 
changes sign, when the screening length in the ionic solvent exceeds a critical value. This leads to 
condensation when the chemical potential is negative. In this state, it is shown using the mean-field 
approximation that spherical and toroidal condensed shapes can be obtained. The thickness of the 
toroidal polyelectrolyte is studied as a function of the parameters of the model, such as the ionic 
screening length. The predictions of this theory should be amenable to experimental verification. 



I. INTRODUCTION 



A distinguishing feature of biopolymers is the presence of charges along the chain and their 
subsequent interaction with each other and ionic aqueous solvents. Yet another aspect is the 
semi-rigidity of some biopolymers. In this sense, biopolymers differentiate themselves from 
polymer melts, which are neutral, and most theoretical treatments involve the treatment of 
flexible polymers. Furthermore, these theoretical treatments of melts almost exclusively use 
models of a self-excluded volume nature. A general treatment of electrostatic interactions 
between charges on a polymer is crucial in order to extend the validity of theoretical models 
to bio-polymers, polymers exhibiting counter-ion condensation in ionic solvents, etc^^iM. 

This paper lays a field-theoretic foundation for a theoretical treatment of charged poly- 
mers in ionic solutions. Technically, the theory applies to flexible polyelectrolytes. However, 
we have argue that by choosing the Kuhn length to be the correlation length of a semi-rigid 
polymer, it should be possible to address the physics of semi-flexible polymers. The intent 
is not only to understand the behavior of the radius of gyration with chain length, but also 
to gain a basic knowledge of condensation phenomena. Our calculations of the radius of 
gyration generalize standard treatments for melts. A mean-field approach was employed to 
develop a phase diagram for condensation/non-condensation of charged polymers. In this 
endeavor, we found the recent experimental paper by Butler et al^ very useful. 

For the case of excluded volume interactions in polymer melts, it was shown by deGennes^ 
and othersiSiiiii^ that in the limit of low number concentration of monomers, the physics of 
polymers is analogous to the onset of a second order phase transition. It will be shown that 
a similar situation ensues when considering polyelectrolytes. The Flory exponent in this 
case, for long chains, is argued to be identical to the classical case. This is in distinction to 
the arguments^ given by Katchalskyi and FloryS. for polyelectrolytes, where a long-range, 
unscreened Coulomb interaction was employed and the effects of any solvent were ignored. 

This paper develops a model which begins with a screened Coulomb interaction between 
native electron charges and adsorbed counter-ions along the segments of a Gaussian chain. 
This screening length is independent of the solvent properties. These charges also interact 
with the surrounding ionic solvent, treated as a continuum. The theory in this paper contin- 
ues the development of a functional integral technique^^ which begins with the formulation 
of a system of many homopolymeric flexible chains, and whose segments interact through a 



finite-range potential. 

The functional integral technique utilized in this paper is different from the coarse- 
graining method employed by Fredrickson et al.J^ Fredrickson et al use the number density 
of the polymer segments as the order parameter. The advantage afforded by this technique 
is that the interaction term becomes quadratic and hence trivial to calculate. However the 
entropy, or free-chain part of the energy functional becomes non-linear, complex and hence 
presents a computational challenge. The order parameter introduced in our paper is a prob- 
ability amplitude whose absolute square corresponds to the number density of segments. 
The free-chain portion of the energy functional is linear in this approach, while the inter- 
action term remains understandably non-linear and difficult computationally. The positive 
semi-definite form of our energy functional may also offer some computational advantages. 
Our formalism can be viewed as a generalization of the self-consistent field theory used by 
Curro's group.^^^ The model of interaction used in this paper differs somewhat from that 
utilized by Curro et al. Our approach is most similar to that of Kleinerl^^J^, but the scope 
of our study is broader. 

Functional integration is used to integrate over the fluctuations in the ionic solvent, 
and derive an effective interaction between the charges on the chain. This interaction is 
repulsive at short distances, and attractive for larger separations. The chemical potential 
/i corresponding to this interaction can be computed. The chemical potential changes sign, 
depending on the parameter regime. Using the interpretation of /i as the energy required to 
add a segment (monomer) to the system, one associates fi > with melt-like behavior, and 
/i < with a condensed state. Indeed it will be shown that when /i < 0, it is possible for a 
single strand of polyelectrolyte to condense into a toroidal shape. Spherical shapes are also 
possible. 

II. THE BASIC FORMALISM 

We reproduce here briefly for completeness, the functional integral formulation of a ho- 
mopolymer chain. For the record, it must be stated that this formalism differs from the 
standard mapping of the self-avoiding walk of a single chain onto a (f)'^ field theory.i^^iS It can 
be generalized in a straightforward way to deal with many chains. Technically speaking, the 
model is restricted to flexible chains. However, if one treats the segment length as a Kuhn 



length, i.e. as an effective distance which represents the correlation length of a semi-flexible 
polymer, then one expects this model to reproduce the overall features of a semi-rigid chain 
as well. 

The probability distribution 60(1, 2; n) for a single Gaussian chain may be represented 
byii: 



60(1, 2; 




where b is the Kuhn length of the polymer as discussed above, and where (9„' = This 
expression is obtained by considering only the entropy of a flexible chain. 

Go{l,2;n) can also be thought of as the Green's function for the diffusion operator in 
three dimensions, as indicated by the last of Eqns. H 

Alternatively, one knows from methods in functional integration that—: 

/3J-'o = J dn'd^x 'ijj*{x,n') 

where V^ip = Vip*Vip, [3 = kB is Boltzmann's constant and T is the temperature. 
Thus we have another way of thinking about a system of flexible polymers, in terms of a 
probability amplitude ip{x, n) and an energy functional jSJ-'o which is isomorphic to one that 
describes diffusion. Here (x, n) labels the location x in physical space, of the n-th segment of 
a chain. A similar approach has been advocated by Kleinert.— Following the convention 
in quantum field theory, n)p is interpreted as the probability of finding a polymer 
segment at a given location in space. 

The main advantage of the functional path integral formalism is that one can model more 
easily interactions in systems with large numbers of polymers by adding an interaction term 
as shown below: 
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f3AJ^ ^ J d^xcfx' dndn'{V{ilj{x,n),ilj*{x,n')) (3) 

where V{ip{x,n),ip*{x' ,n')) represents the interaction between segments. Note that in gen- 
eral, on physical grounds, n ^ n', as different configurations of the strand permit segments 
from different parts of the chain to be in the proximity of each other. 
The form chosen for V in this paper is: 

V(^(a;,n),V^*(f,n')) = Q mx,n)\' Ur,,n' {\x - ^ \) ,n')r (4) 

Here Un,n'{x — x') represents a repulsive-attractive short-ranged interaction between 
monomers. The precise form of the potential useful for describing polyelectrolytes will 
be discussed in detail the next section. 

Upon extremizing the function j3T with respect to ip* , one obtains a non- linear diffusion 
equation: 



V^il){x,n) + J dn' J d^x'\ip{x' ,n')\^Un,n'{\x - x'\) '4j{x,n) = (5) 

This equation is also analogous to the single-electron Hartree-Fock equation in quantum 
mechanics, and applies to a single polymer strand. 

If a generalization to many strands (several electron orbitals in quantum mechanics) is 
required, then the requisite form would be: 

-jV^-^i{x,n)+J2 I dn' J d''x'\^j{x\n')\^Un,n'{\x-x'\)^,{x,n) = (6) 

where the subscript on the amplitude indicates a chain label, Nc being the total number 
of chains in the system. The rest of the paper will focus primarily on the physics of a single 
chain. 

We shall now show how a connection can be made to standard results in polymer theory. 
The correlation function Qo{k,p), p being the separation between any two segments on the 
Gaussian chain can be obtained by using the Fourier transform of the Gaussian distribution 
function discussed earher in the section. This Fourier transform shall be denoted by (?o(/c, a;): 
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A simple contour integration, performed by closing the contour in the complex uj plane 
in the clockwise direction leads to: 

go{Kp)^eM-\p\ {bkf/6) (8) 
The radius of gyration -Rq(p) can be defined as: 



/ k=0 



— — power of k = 2 (9) 

The experimentally accessible structure factor associated with the Gaussian chain of N 
segments can defined in the following fashion: 



So{k) = N-^ / dmdn Qoik, \m - n\) 
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^ N {1- (1/3) k''Rl{Nf) + 0{k^) (10) 
Rg{p) can also be defined in terms of the structure factor: 
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— = power of k = 2 (11) 

In the presence of interactions, Go(^,<^) gets renormalized in the usual way, and an 
associated dressed correlation function Q{k,N) and dressed structure factor S{k) can be 
defined analogously: 



g{k,N) = [ ^exp{-tujN) G{k,uj) 

S(k) = N-^ g(0,N)-^ / / dm dn g(k,\m-n\) (12) 
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where a normalization factor has been included in the definition of the structure factor. 

The radius of gyration Rg{N) in the presence of interactions can be defined following the 
ideas for the free chain in terms of the dressed correlation function: 
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where a normalizing factor has been inserted on the right hand side of Eqn]13lfor complete- 
ness. Equivalently: 
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III. POLYELECTROLYTE IN A SOLVENT 



Let us suppose that the net charge density on any segment along the chain is given by 
sum of negative and positive charge densities on the chain: pt{x,n) = p^{x,n) + p^{x,n), 
with p_ < 0, p+ > 0. Then the total free energy which describes the interaction between 
these charges and the electrolytic solution is given by three terms, describing the interaction 
of charges along the chain, the solvent, and the interaction between the polymer and the 
solvent, respectively: 



P 3 ~l~ polymer ~l~ solvent ~l~ polymer— solvent^ 




(15) 



where 0(x) is the electrostatic potential of the ionic solvent, Kg = l/Ag is the inverse screening 
length of the polymer, k = 1/A is the inverse screening length in the solvent, where e is the 
dielectric constant of the solvent (e ~ 80 for water). The reason for assuming the screening 
lengths to be distinct for the chain (Ae) and the solvent (A) is that even in the absence of an 
ionic solvent, when the adsorbed charge is null, the conformations of the negatively charged 
Gaussian chain can screen the Coulomb interaction between the native charges on the chain. 
The charge distribution on the chain is given by: 



Here q-{n) refers to the native negative charge on the backbone of the polymer, and 
g+(n) refers to the positive counter- ions that may be adsorbed onto the chain. One can 
conclude from this that if the negative native charges are singly charged, then the adsorbed 
counter-ions would have to be multiply charged, as otherwise the net charge would be zero, 
and we would be reduced to studying a neutral melt. It also allows for the possibility that 
the counter-ions are docked in a physically different location along the chain than the native 
charges, which would allow for local regions of attraction to be created. 

Beyond these inferences, one notes that if the above function Eqnll5lwere to be extremized 
with respect to the the scalar potential 0, one obtains the standard Debye-Huckel equation, 
valid in the limit of low ionic strengths. Thus our paper goes beyond the standard, mean- 
field Debye-Huckel model by considering fluctuations around that approximation in the full 
energy functional given by Eqn. ^1 These fluctuations will be shown shortly to lead to an 
attractive component of the correlations between like charges along the chain. 

Note that the electrostatic potential appears quadratically in the functional. One can 
now integrate over solvent degrees of freedom, in particular, accounting for deviations from 
the Debye-Huckel mean-field limit to obtain an effective functional: 
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FIG. 1: This is a plot of the effective potential derived in the paper. Representative parameters 
were used: e = 80, Ae = 3AA, A = lOl. 
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One sees that the effective potential V, exemplified in Fig JIIII is given as follows: 



V{\x-x'\)= ^^v(-^e\S-x'\) _exv{-.\x-x'\) 
\x — x'\ e \x — x'\ 

Observe that integrating over the solvent degrees of freedom has led to an attractive (neg- 
ative) supplement to the original short-ranged screened Coulomb interaction. The electrical 
polarization of the water, signified by e serves to weaken this attractive part of the potential. 
Given that the solvent we consider is water, which has a natural pH of 7, one always has an 
associated non-zero screening length viz., A. 

The various parameters used in model, e.g. q'+(?T.), A, Ae are inter-related. This is because 
changing the ionic concentration in the solvent affects not only A, the solvent screening 
length, but also the amount of adsorbed charge q+{n) and hence the screening length Ae 
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FIG. 2: A schematic view of the correlation effect which leads to an effective attraction between like 
charges on the chain. Opposite charges in the solvent create a region which attracts like charges. 

along the chain as well. And it may be prudent initially to appeal to experiment in order 
to assess their magnitudes. 

The physical notion behind the effective potential is analogous to what happens in 
electron-phonon physics in metals, when the attraction between electrons and phonons (mo- 
tion of positive ions) leads to an effective attraction between the electrons themselves. This 
effect can be depicted pictorially in Fig lllll 

The approach taken in this paper, viz., integrating over the solvent degrees of freedom, is 
complementary to the approach of Manningi^, in which the focus rests on the electrostatic 
potential of the solvent. The polyelectrolyte, represented by an infinite line charge, was 
shown by Manning to polarize the surrounding ionic solvent. There exists a thin annular 
cylinder around the rod which attracts counter-ions, allowing for condensation under certain 
conditions. Manning's theory is a mean-field approach. The approach developed in this 
paper is more general. The present theory accounts for fluctuations around the mean-field 
approximation. It allows for a general description of the physics of polyelectrolytes, and is 
not restricted to a description of condensation alone. 

Furthermore, as suggested earlier, if one takes the segment length b to represent the 
correlation length of a semi-rigid polymer, we expect the model to reproduce the overall 
phenomena associated with semi-rigid polyelectrolytes. 



IV. REGIMES OF POLYELECTROLYTE BEHAVIOR 



The effective potential derived in the previous section permits an insight into the variety 
of behavior that polyelectrolytes in solution can display. In order to do that, it is useful 
to add a constraint to the energy functional described in the previous section, viz., Tefj to 
conserve the number of monomers on the polymer: 

/3AJ^^ =-^1 j dn j d^x\ilj{x,n)\^ 
/i = 27rA^4co(A^ - A^e"^) 

4 =P{q+-\q-\f (19) 

where £b is the Bjerrum length, cq is the number density of monomers, is the chain 
length, and the chemical potential is obtained by minimizing the energy functional in EgnlTTH 
in the limit that the wave amplitude if) is independent of x, n. It has also been assumed for 
simplicity that there is a uniform charge distribution along the chain. 

By convention, the chemical potential ^ represents the energy required to add a segment 
to the system. Hence, upon examining Eqn ll9t one sees three possible regimes. 

Melt: When, /i > 0, the polymer behaves as a melt, when repulsion dominates, and it 
costs energy to add a segment to the system. This occurs when the screening length along 
the chain Ae > A/y^. This is similar to Manning's discussion of polyelectrolytes. 

O point: The monomers are perfectly miscible with each other when the chemical poten- 
tial /i = 0. This occurs when Ae = A/ 

Condensation: When the screening length along the chain Ae < A/y^, the chemical 
potential becomes negative, indicating that the system encourages the addition of monomers. 
It is natural to identify this regime as that of condensation. Our criterion is similar to 
that of Manning^^, but it requires a consideration of fluctuations around the Debye-Huckel 
approximation. In this sense our approach bears a greater resemblance to the analytic 
arguments in the paper of Jensen et al^. 

As discussed earlier in the paper, if we assume that the native charge on the chain is 
negative, viz., g_ < and corresponds to the charge of a single electron, then clearly, we 
require g+ the charge on the condensed counter- ion to be multi-valent, as otherwise the net 



charge is null, leading to a neutral polymer, and a discussion of condensation becomes moot. 
It is beyond the capability of the current theory to address the difference in the condensation 
behavior of different types of multi-valent molecules carrying the same charge^. 

Using SAX techniques, Butler et al^ show that virus solutions need a minimum diamine 
concentration for the onset of condensation, as signaled by the emergence of a peak at a non- 
zero wave-vector. Moreover, they show that there also exists a maximum ion concentration 
above which condensation ceases. Equation^Jcan be interpreted in light of these experimen- 
tal results. For low ion concentrations, g+ the adsorbed counter-ion charge is correspondingly 
small, and the screening length along the chain Ag is very large. Since the dielectric constant 
of the aqueous solvent is high (~ 80), the chemical potential fi = 27rA^4co(Ag — e^^ A^) can be 
greater than zero, preventing condensation. As the ion concentration increases, the screen- 
ing lengths along the chain and in the solvent decrease. When this change occurs in such 
a way that Ae shrinks faster than A, such that A > eA, then the chemical potential yU < 0, 
leading to condensation. Indeed this criterion is similar to the one derived on experimental 
grounds by Butler et al. For even higher ionic concentrations, Ag cannot shrink much below 
the monomer length, and it is possible for the chemical potential to become positive once 
again. Condensation then ceases. It may prove possible to adjust the parameters in this 
theory to reproduce quantitatively the condensation data of Butler et al^. Our aim here 
is to display that the theory is simply capable of addressing experiments. The cessation of 
condensation at relatively high ionic concentration appears to be a new experimental result, 
and it is important to point out that the theory in this paper can address this issue. 

The theory presented in this paper is similar to that of Golestanian and Liverpool^i. 
These authors use a phenomenological theory of semi-rigid chains decorated with charges. 
The only interaction in their theory is that between these charges. The presence of a solvent 
is implicitly acknowledged through an effective screening length. In our theory, it would 
correspond to neglecting the terms T solvent and T polymer -solvent- They show that fluctuations 
in the shape of a semi-rigid rod lead to a buckling instability. These authors too find 
cessation of condensation upon increasing the salt concentration sufficiently, by considering 
the behavior of the correlation length. 

The current theory can accommodate the zipper-like motif invoked by Kornyshev and 
Leikeii^^. Our approach can be construed as providing an analytic background to the nu- 
merical simulations of Stevena^^ of bead and spring polyelectrolytes in a solvent. 



The onset of condensation in biopolymers, based on the assumption of entirely rigid 
polyelectrolytes, has been shown by Kholodenko^^ to be analogous to the melting of vortices 
in Thouless' two-dimensional model of a glass. In the future, one could attempt to investigate 
aspect further, using the current formalism, which does not assume rigid rods, as follows. 
Rather than integrating over the the electrostatic potential of the solvent, one could attempt 
a functional integration over the polymer degrees of freedom. In this case, extremizing the 
new effective functional with respect to the electrostatic potential will lead to a generalization 
of the Kholodenko model to non-rigid polyelectrolytes in solution. 

Co 0+: In addition to the three regimes discussed above, there is one more important 
regime within the melt phase which merits attention, and permits a connection with the 
scaling arguments first derived by de Gennes^. Consider the energy density Pfeff for a 
homogeneous ip: 

Pfeff = il/2)NaV{k = 0)m^ - col^n (20) 

where a = ib/b, where the dimensionless number concentration cq can be expressed in terms 
of the actual monomer number concentration Cq and the Kuhn length 6 as Cq = Co(&/-\/6)^ 
Thus when co(&/-\/6)^ << 1, one sees that the double- well structure displayed in the next 
figure is just beginning to form. This is the traditional signal for a second-order phase 
transition (see Fig llVjl . Physically this occurs because the mean free path between collisions 
is very large in the low number concentration limit, permitting large fluctuations, the hall- 
mark of second order phase transitions. Note that the criterion above encompasses the 
arguments given by Freed^'', where it is shown that phase transition-like behavior occurs in 
the 6 — > limit. 

V. RENORMALIZATION GROUP CALCULATIONS IN THE MELT REGIME 

We shall now show that the polymer model considered in this paper exhibits a fixed point 
of the Renormalization Group transformation, reflecting the large fluctuations expected in 
the cob^ 0+ (melt) regime. Since a Renormalization Group fixed point in the long wave- 
length limit (for long chain lengths) is always independent of the details of the interaction, 
it is expected that our model of realistic interactions should yield results identical to the 
excluded volume model usually considered in polymer theory. Nevertheless, it is useful to 
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FIG. 3: A schematic of the double well potential given by the energy density Pfeff, in arbitrary 
units. Notice that the height of the potential barrier decreases as the concentration of monomers 
decreases, allowing for more fluctuations between the minima. 

derive this result explicitly so that it displays the validity of the field theoretic approach 
employed in this paper. Further, we shall compute a reasonably accurate value for the Flory 
exponent, obtained when the chain length is very large. In what follows, we shall use a 
natural system of units, in which the unit of length will be taken to be b/\^, where b is 
the segment length. As discussed earlier in the paper, b is the monomer length for a flexible 
chain, b could be taken to be the correlation length for a semi-rigid chain. 

Finally, note that the model is analogous to the theory of dynamical second order phase 
transitions^^, with the segment label n playing the role of time. The model has four inde- 
pendent variables, viz., three spatial dimensions and an additional label for the location of 
a segment along the chain. 

Towards that end, let us consider the first order correction to the bare vertex (the basic 
interaction displayed in Eqn E)) • It corresponds to the first order polarization diagram in 
many-body physics. It represents the generalization of the correction that occurs in an 
excluded volume interaction model^^. 

The result depicted in this figure can be encapsulated as the renormalization of the 
coupling constant a (using a length scale of b/\/6): 
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FIG. 4: The Feynman diagram indicating schematically the bare vertex, whcih represents the 
interction term in the energy functional. It displays two segments interacting with each other at a 
distance. 
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where the Fourier transform Go{p,uj) of the unperturbed propagator (Green's function) 
defined in Eqn. ^ is given by: 
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Now no(q', uj) represents screening of the interaction between monomers due to intervening 
segments. Thus it is tacitly assumed that we are in the long-wavelength regime, and the 
chain length is large, as it is in this regime that one expects screening effects to dominate. 
Screening is clearly unimportant for extremely short chains, as the probability of finding 
intervening segments in this case is vanishingly small. 

Upon performing the frequency integration in Eqn. |^ by the method of residues, by 
closing the requisite contour in the lower half-plane: 



no(g,w) = - 
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where Tie implies the real part of the quantity on its right hand side is to be considered. 
The polarization IlQ{q,Lj) can be recast as: 



no(g,t^) = - {-^^j ^cos(t^A) J dpp{exp{-Cp) - exp{-2pqX)) 



(24) 



where C ^ 0"*", and the factor of exp{—(p) has been inserted into Eqn. |^ to achieve 
convergence. We now invoke the integral representation of a gamma function, viz., 
/q°° (ixx^ exp(— ax*") = m~^a~^'P~^^^^"^T((p + l)/m) to obtain closed form expressions for the 
integrals in Eqn. |211 However this gives rise to divergent factors such as r(— 2). These di- 
vergences can be cured via different methods, such as Wilson's renormalization scheme. The 
renormalization scheme to handle ultra-violet divergences utilized in this paper follows more 
closely the method of dimensional regularizatioit^P of t'Hooft and Veltman. This technique 
allows one to separate out the infinities, by setting p ^ p ~ e, e 0+, for p + 1/m equal 
to a negative integer, in the expression above for the integral representation of the gamma 
function.. This allows us to express the polarization function as: 

Using the Macclaurin series for r(— 2 + e): 

r(-2 + e) ^ (1/2) {e-' + ^(3) + 0{e)) , (26) 

one can isolate the divergent portions of the polarization function as being proportional to 
and e~^. These infinities can be canceled by the invocation of appropriate counter-terms 
in JF, as is normally done in field theory. This allows us to retain only the finite portion of 
the polarization function: 



Uo{q,u) = - 



/H =-[-^) (ln(kl)-3/2 + 7) (27) 

where 7 is Euler's constant, numerically close to 0.5772. 

A test of the renormalization scheme used to handle ultra-violet divergences here is to 
calculate an observable or a known quantity such as the Flory exponent. This will be done 
shortly, after a thermodynamic fixed point has been identified. This fixed point is associated 
with the second order phase transition discussed in the previous. 

One must go to the long wavelength limit in order to invoke the Renormalization Group. 
This Renormalization Group applies to the infra-red limit (g 0), and must be distin- 
guished from the preceding discussion regarding the behavior in the g ^ 00 limit. 



It is convenient to redefine a scaled coupling constant in the g ^ limit, with an associ- 
ated /5-function which yields the Renormalization Group flow: 



. f f{u;)V{q = Oy 
aR{q,uj) =aji \ 
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where q = L^^. In analogy with the theory of second order phase transitions, scale invariance 
of the coupling constant must be invoked for a proper description of the physics. It follows 
from Eqn. EHlthat P{q:ji) has a non-trivial flxed point at aji = a*j^, or at aji = a*^: 
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f{uj)V{q = 0) 

It follows from Eqn]29lthat in the long wavelength limit, in the Cq 0^ regime, the inter- 
action has the following universal form, independent of the particular interaction potential 
we started off with: 

aViq 0) ^ - \ (30) 

^ [HM) -3/2 + 7) 

The same result could have been obtained in a more mundane fashion by using the so- 
called Random Phase Approximation (RPA) in which one traditionally sums up bubble 
diagrams of all orders. Nonetheless, the Renormalization Group approach is essential in 
that it explicitly manifests the scale invariant properties of the system. 

Scaling arguments will be utilized to obtain a value of the Flory exponent. These are 
provided in lieu of the conventional Renormalization Group argument based on scaling trans- 
formations. That treatment is based on the scaling properties of the simpler 0^ Ginzburg- 
Landau functional. In the present paper, we have a more realistic interaction potential 
viz., V to deal with. Consequently, the straighforward scaling properties of the 0^ excluded 
volume theory are no longer valid. A more complex set of scalings, involving the screening 
lengths are required to maintain the scale invariance of the current theory. But this implies 
a restriction on the class of polyelectrolytic solutions that can be considered, conforming 



to the scaling transformations. The alternate derivation offered below suffices to yield the 
Flory exponent fairly accurately. 

Consider scaling behavior of the energy functional for a Gaussian chain: 

P^o = J^ j ^,r{Kuj){-iu: + k')^{k,u:) (31) 

The diffusion-like, self-similar form permits Functional to be scale invariant under the 
following transformations: 



k ^k/i{M) 

^ ^^^5/2(^) (32) 

where uq = 1/2. The relation between ^ and A/" is identical to the relation between the 
radius of gyration Rg and A^. Hence we shall identify the length scale ^ (A/") with the radius 
of gyration. 

Let us now consider the case when the interaction viz., V is turned on between the seg- 
ments. The behavior of the energy functional in the long wavelength limit {q 0,uj 0) 
is dominated by the behavior of the effective interaction at the Fixed point of the Renor- 
malization Group transformations (see EgnlHUj): 

^•-^/^/S^-f*--'!^);'*'^' '''' 

where we have used x'^ ~ 1 + elnx, e << 1 to convert a logarithm to a power law. 
This functional exhibits invariance under: 



k ^k/i{M) 

tP -^4^e^\Ar) (34) 



where iy = 2/3- l/{3n^) ^ 0.63. 



FIG. 5: A diagrammatic depiction of the approximation used in the calculation. The thick line 
indicates G. The first term on the right hand side of the equation refers to Go- The second term 
refers to the exchange-like diagram referred to in the text. The final term is the second order 
correction. 

This estimate of the radius of gyration holds for asymptotically large chain lengths A^, 
when the concentration of monomers Cq << {V6/bY. 

If one accepts the value of z/ ^ 0.5889 as obtained within the framework of the Edwards' 
model, then the result in this paper is accurate to about 92%. On the other hand, H. 
Kleinerip^ has shown that when the asymptotic series in the e expansion is handled appro- 
priately, the value of 0.5889 is replaced by 3/5. Kleinert then uses field theoretic methods 
and the replica trick to rederive a value of z/ ~ 0.62. The value we obtained, u ~ 0.63 is 
fairly close to the universal value for u ^ 0.634 accepted in phase transition theory. 

This value is quite different from the one derived by Flory® and Katchalsky^ for the case of 
an unscreened Coulomb interaction, which yielded a linear dependence on the chain length. 

VI. RADIUS OF GYRATION FOR SHORT CHAINS 

For the case of short chains, when the coupling constant a = ^/6ib/b < 1, one can 
compute the lowest order Feynman diagram to estimate the self-energy, which renormalizes 
the Green's function GQ{k,u!). This estimate is correct when a < 1 and when the chains 
length is relatively short. This is because screening is expected on physical grounds to be 
unimportant for short chains. It turns out that there are two terms of 0{a) in a perturbative 
expansion of the self-energy. The first one, the tadpole diagram can be renormalized away 
in the usual fashion by the addition of a counter-term. The remaining term is reminiscent 
of the exchange diagram in many-electron physics (as shown in Figur^. 

Using the Feynman-Dyson technique, repetitions of this diagram can be summed as a 
geometric series to infinite order. This exchange diagram contribution to the self-energy can 
be expressed in closed form as follows: 



J^^'^oi^^ - q\) (35) 

where it is tacitly assumed that we are in the low concentration limit. Since our focus is 
on estimating the radius of gyration for relatively short chains, we need only to evaluate 
the self-energy in the long-wavelength limit, viz., k 0. To facilitate this, the angular 
integrals in EqnlH^l can be done and the contributions to 0{k'^) of the integrand can be 
written explicitly: 
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(36) 



V 3 ; [e^q^ + K^y (g2 + /t,2) 

The subsequent fc— integration can be performed with the aid of Mathematica. The results 
are long and cumbersome and not much is gained by stating the expressions explicitly, other 
than to say that the integrals contain logarithmic terms of the form Intu. Formally, the next 
step is to evaluate the correlation function: 



exp(— icuA/") 



— (37) 

-iuj + P - T.exch{k,u;) 

This integral could be evaluated using the method of residues, by closing the contour in 
the lower half plane if the root(s) of the denominator could be located. Since our interest is 
in the k ^ Q limit, it is possible to estimate the root of the denominator perturbatively: 



^root 



-ik^ + iT,exch{k, uj = —ik"^) + O(a^) 



(38) 



Before the residue can be calculated, the presence of the logarithmic terms referred to 
above imply the existence of a branch cut in the complex uj plane from — oo to along 
the real axis. Hence the contour has to distorted slightly into the lower half-plane along 
the negative real axis to avoid this branch cut, so that the only singularity enclosed by 
the contour in the lower half-plane is at the root given by Eqn]38l It then follows, using 
Mathematica, that in the long wavelength limit: 



Radius of gyration ~ N 



FIG. 6: Schematic depiction of the hnear behavior of the radius of gyration for short chains. Note 
that the Kuhn length can be shorter than the segment lenght due to a fluctuation effect. 



g{k, N)/g{k = 0,N) ^ exp{-ikg{N)) 

g{N) =a{Xt- A^)r^ (39) 

where Ae and A are the screening lengths along the chain and in the solvent, respectively. 
The experimentally measurable structure factor can then be evaluated: 



S{k) ^ N{1 - {l/3)eRl{N) + 0{k^)) 

= ineffective /VQ) N 

-- ilSa) 



^ef fective 



(40) 



Note that Rg{N) ~ in this regime, as displayed in Fig lVIl Moreover, it turns out 
that Rg{N) is extremely sensitive to the parameters in the potential. This is in contrast to 
the asymptotic regime encountered in the previous section, where universal behavior was 
found. X-ray or neutron scattering experiments could be performed to test the predictions 
presented in this section. 

One could interpret beffecUve in EgnHUl as an effective correlation length and suggest 
^effective ^ ^ aS a SI gnal for buckling, and hence for condensation. This interpretation is 
equivalent to the chemical potential considerations presented earlier in the paper, and is 
similar to the approach of Golestanian and Liverpool^i. 



VII. FINITE CONCENTRATION OF SEGMENTS 



The discussion in the preceding two sections has been primarily focused on behavior in 
the melt regime, when the number concentration of segments is vanishingly small. In this 
case, fluctuations were accounted for, and the radius of gyration for both long and short 
chains was computed. As the concentration of segments increases, the segments become 
packed increasingly closer, thereby decreasing fluctuations in the system. Mean fleld ap- 
proximations, obtained by extremizing the energy functional can then be invoked to obtain 
insight into the physics. 

In a previous paper,— which utilized an excluded volume interaction model, we showed 
how tube-like structures can be obtained. If one restricts attention to obtaining an enve- 
lope of structures obtained in the various regimes delineated in section III, then a similar 
technique provides useful insight in the current approach as well. The basic idea, designed 
to ease computation, is to replace the short-ranged potential by an effective delta-function 
pseudo-potential. The effective coupling constant which characterizes the pseudo-potential 
can be positive or negative, depending on the value of the chemical potential /i. As discussed 
earlier, the value of the chemical potential is an average way of determining whether the 
attractive part or the repulsive part dominates the behavior of the system. 

The advantage of this method is that it yields the correct average behavior of the system 
for a reasonably small effort. The disadvantage is that if one is interested in details of the 
structures which change on the spatial scale less than the one over which the interaction po- 
tential varies, then one must resort to vastly more detailed calculations. The approximation 
consists of the following replacement: 




(41) 



Extremization of the functional leads to: 



(— - + aV(k = 0)\'^(x,n')\^ - fiU(x,n) = (42) 
on 

For cases when the segment label n is not physically relevant, this equation reduces to: 



(-V^i (|^(x)p-l))^(a;) = (43) 
where the distances are scaled in units of (6/ VQ){a\UMorse{k = 0)|)^/^, and the amphtude 

— 1/3 

ip has been scaled by Cq so as to be dimensionless. The positive sign in the non-linear 
partial differential equation holds when yU > in the melt phase, and the negative sign when 
in the condensed state (/x < 0). 

Associating yU < with a condensed state implies the existence of coherent structures. 
The first structure we investigate is one in spherical geometry, with a larger amplitude at the 
center than towards the edges. The ordinary differential equation that requires a solution 
is: 

-;^|:('-'^)-(l*('-)l^-iW'M = o (44) 

It turns convenient for numerical purposes to use y = 1/r, so that EqnJ44l becomes: 

-y^^"{y)-{^{yf-l)^{y)=0 (45) 

For y 0, which is the same as r — > oo, the wave amplitude is expected to decay, so 
that the non-linear term vanishes, and yields "^{y) = ay exp{l/y). One can now integrate 
numerically from some y = y minimum to y = y maximum = r minimum, adjusting a such that the 
slope of the amplitude vanishes at y maximum- The result is shown in Figure IVlTl 

Another interesting structure that we have investigated is a toroidal structure. The 
interest in this structure arises from the fact that strands of DNA in an ionic solution 
condense into such shapes under appropriate solvent conditions^LSS. One may construe the 
model above (Eqn HH|) as an effective model for DNA in solution. The co-ordinate surfaces 
in toroidal geometry are^°: (i) planes through the z-axis, represented by an azimuthal angle, 
(ii) spheres of varying radii centered up and down the z-axis, and finally, (iii) tores, or anchor 
rings around the z-axis, labeled by the location of their centers at a distance a, and cross- 
sectional and axial radii a cschr/ and a cothr/ respectively, for r] = constant, rj being the 
toroidal co-ordinate. Since we seek solutions with toroidal symmetry, the only independent 
variable we need to consider is rj, which leads to a non-linear ordinary differential equation: 



2 smh{r]/2f{tanh{T]/2)ij'{r]) - 2 sinh(r//2) - a^ii^ivY - l)V^(^) = (46) 
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FIG. 7: Displayed is the radial profile of a spherical coherent structure obtained for fi < 0. The 
co-ordinates are dimensionless. 

The radius a of the torus is determined self-consistently as the value which yields a zero 
slope but non-zero amplitude at rj ^ oo. This is the center of the doughnut. Again, since 
we are looking for a doughnut shaped object, it follows that the slope of the function should 
be zero at a value of rj at which the cross-sectional radius of the tore is a cschrj = a, i.e., 
when 1] = arcsinhl ^ 0.8813. At this point the the z-axis is a tangent to the tore at the 
origin, and physical consideration implies a zero slope for continuity. The value of the wave 
amplitude is taken to be zero at rj = arcsinhl, i.e., the origin. A shooting method was 
employed where a was varied iteratively until a solution with a zero slope at ?7 ^ oo was 
obtained numerically. Operationally, the equation was integrated to some large value of t]. 
The result is displayed as contours in Figure VIII in the x — z plane.. In effect a doughnut 
shaped structure is obtained, whose hole is partially filled. The energy of this structure is 
identical to that of the spherical blob displayed in Fig. I VI II 

These are two examples of coherent structures that are possible for fi < 0. They are 
both equally energetically favorable. As such, the approximations employed during these 
calculations are applicable for large chain lengths, when it is known exoperimentally that 
either spheroidal or toroidal configurtions are equally likely to occu r^^i^^ . These examples 
do not constitute an exhaustive list. 

Finally, we note from the previous figure, displayed in dimensionless variables, that the 
width of the toroidal configuration is 0{1), in the length scale used. Based on this observa- 
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FIG. 8: Displayed are the contours of a toroidal structure for /x < 0. Note that the doughnut has 
a profile which indicates a changing density with distance from the center. Note that the width of 
the structure is approximately unity in dimensionless units. 

tion, one can calculate the thickness of the toroid for various parameters, and an example 
is given in Fig. IX. Experiments of the sort performed by Golan et al2& and by Butler et al^ 
should be able to provide experimental verification of our predictions. 



VIII. CONCLUSION 



The theory described in this paper develops a functional integral technique to treat re- 
alistic interactions between segments of a polymer in a realistic way, through the use of 
a finite-ranged repulsive-attractive interaction potential. Examination of the chemical po- 
tential led to a classification of homopolymeric systems. It was pointed out that such 
a classification would be impossible with an a priori excluded volume interaction model. 
Renormalization Group techniques were used to show that standard concepts in polymer 
physics are recovered in the limit of low monomer number concentration, for asymptotically 
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FIG. 9: The width of the toroid is plotted in nano-meters for two different cases, as the chain 
length is varied, for two different screening lengths in the solvent. The upper curve coresponds to 

a screening lenght of 50b/sqrt6. The lower curve coresponds to a screening lenght of 10b/sqrt6. 

long chains, in the melt state. The radius of gyration for extremely short chains is also 
calculated, and is linear in the chain length, reminiscent of a semi-flexible chain. When the 
chemical potential is negative, condensed structures are shown to exist, both in spherical 
as well as toroidal geometry. The predictions that follow from the theory presented in this 
paper, viz., the short chain radius of gyration, the widths of the toroidal configurations as 
functions of various experimentally accessible parameters could be verified experimentally. 
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